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ABSTRACT 

Context. Supernovae are known to be the dominant energy source for driving turbulence in the interstellar medium. Yet, their effect 
on magnetic field amplification in spiral galaxies is still poorly understood. Previous analytical models, based on the evolution of 
isolated, non-interacting supernova remnants, predicted a dominant vertical pumping that would render dynamo action improbable. 
Aims. In the present work, we address the issue of vertical transport, which is thought to be the key process that inhibits dynamo 
action in the galactic context. We aim to demonstrate that supernova driving is a powerful mechanism to amplify galactic magnetic 
fields. 

Methods. We conduct direct numerical simulations in the framework of resistive magnetohydrodynamics. Our local box model 
of the interstellar medium comprises optically-thin radiative cooling, an external gravitational potential, and background shear. 
Dynamo coefficients for mean-field models are measured by means of passive test fields. 

Results. Our simulations show that supernova-driven turbulence in conjunction with shear leads to an exponential amplification 
of the mean magnetic field. We found turbulent pumping to be directed inward and approximately balanced by a galactic wind. 

Key words, turbulence - magnetohydrodynamics - ISM: supernova remnants, magnetic fields 



1. Introduction 

Large scale magnetic fields, as observed in disk galaxies, are 
thought to be the result of a turbulent dynamo. Mean-field mod- 
els based on the turbulent c-effect can indeed explain most 
of the observed features in many nearby disk galaxies (cf. 
lBecketalJl996h . However, uncertainties remain in the estima- 
tion of the turbulent transport coefficients and the nonlinear sat- 
uration process. Additionally, the helicity conservation in ideal 
magnetohydrodynamics (MHD) could lead to the so-cal led 
catastrophic quenching (Brandenburg & Subramanian 2005), a 
saturation process that scales with the magnetic Reyn olds num- 
ber. H elicity transport by a wind, as considered by ISur et alj 
d2007h . would be one solution to this scenario. 

Recently, corrobative evidence has been found that the de- 
gree of regularity of m agnetic fields decreases with increasing 
star formation activity (Chyzy 2007). This has renewed the the- 
oretical interest in dynamo models based on the interaction of 
supernova-driven turbulence with galactic differential rotation. 

Despite the success of mean-field models, until now there 
has been no clear numerical verification of the dynamo process 
in the galactic context, and semi-analytical models, which are 
based on doubtful assumption s, only have limited predictive 
power (see lGressel et al.|[2008l for a recent account). 

iBalsara et alJ d2004l) found growing small scale 
fields for a non-rotati ng isothermal gas. Similarly, 
Ide Avillez & BreitschwerdJ d2005l) considered the most 
realistic model of the interstellar medium (ISM) to date, 
but did not include the rotation and shear necessary for a 
mean-field galactic dynamo to operate. Except for t heir neglect 
of th ermal instability (TI), the simulations of iKorpi et aU 
(1999) are very similar to ours. However, due to the limited 
computational resources available at the time, their setup suf- 
fered from a too small box size, which prohibited a long-term 
evolution into developed turbulence. 
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Parker] d 19921) proposed the magnetic buoyancy instability 
supported b y cosmic rays as a driving m echanism for the galac- 
tic dynamo. Hana sz"et alj d2004l 12006}) performed simulations 
of this process, employing the injection of cosmic ray energy 
into an adiabatic gas as the driving mechanism. In these simu- 
lations, cosmic ray transport is treated in the (anisotropic) dif- 
fusion approximation. The authors found an amplification of 
the mean magnetic field with an e-folding time on the order 
of 100 Myr. Contrary to these cosmic ray models we consider 
the injection of thermal energy of clustered supernovae and ap- 
ply a radiative cooling function to reflect the multi-phase na- 
ture of the ISM. Although currently one cannot yet resolve all 
the complex effects of the interstellar plasma, going beyond 
simple artificial forcing (relying on ad-hoc assumptions), these 
new approaches mark the first step towards realistic dynamo 
models based on fundamental physics. 



2. Methods 

We explore the evolution of the galactic magnetic field within 
the stratified, turbulent ISM by means of combined 3 D, resis- 
tive M HD simulations employing the NIRVANA code (Ziegler 
2004) and ID mean-field models based on turbulence parame- 
ters. 



2. 1 . Direct simulations 

The adopted domain covers a vertically centred box of 0.8 x 
0.8 x 4.267 kpc 3 , representing a local patch of the galactic 
disk. The physical model includes an external gravitational 
potential dKuiiken & Gilmora Il989l). differential rotatio n in 
the shearing box approximation dGressel & Zieg ler 2007), and 
optic ally- thin radiative cooling. For the latter, we apply a sim- 
ple piecewise power- law description dSanchez-Salcedo et al.l 
120021 ISlyz et al.ll2005l) . In the following, we will briefly intro- 
duce the key parameters of our simulations. For a more detailed 
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descrip tion of the model, we refer the reader to iGressel et alJ 
(120081) . 

We choose an initial midplane density of «o = 1 cirT 3 at a 
mean molecular weight of fi = 0.6. This sets the equilibrium 
midplane pressure at a value of po/ks = 6000 KcrrT 3 . The ini- 
tial stratification is an exact hydrostatic equilibrium solution in 
the absence of the kinetic pressure caused by the SNe. Contrary 
to previous isothermal solutions, our approach applies an effec- 
tive equation of state given by the additional constraint that the 
initial disk is in equilibrium with respect to the radiative cool- 
ing. The weak initial magnetic field has a radial and azimuthal 
component with B^/Br = -10. The magnitude of the field as 
a function of height is scaled with (p(z)/ po) 1 ^ 2 to yield a con- 
stant plasma parameter of Bp = 2xl0 7 throughout the disk. The 
rotation period is measured in units of Qo = 25kms~ 1 kpc~', 
while differential rotation is quantified via the shear parameter 
q = dln£2/dln/?. We apply a constant kinematic viscosity of 
v = 5x 10 24 cm 2 s~' along with a microscopic magnetic diffu- 
sivity of 77 = 2x 10 24 cm 2 s -1 , resulting in a magnetic Prandtl 
number of Pm = 2.5. For reference, we also define the mag- 
netic Reynolds number Rm = L^Qq/tj 2500. 

Turbulence is driven solely via the localised injection of 
thermal energy. While type-I SNe are exponentially distributed 
with a scale height of 325 pc, for the more frequent type-II SNe, 
we adopt a prescription where the probability distribution is de- 
termined by the vertical gas density profile. This avoids an ar- 
tificial fragmentation of the inner disk observed in conjunction 
with a static SN distribution. The reference explosion frequen- 
cies are cr^ — 4 Myr -1 kpc -2 and cr? = 30 Myr~' kpc 2 , respec- 
tively, and the associated energies are 10 51 and 1.14xl0 51 erg, 
where the higher value for the latter accounts for th e contri- 
bution from a stellar wind of the massive progenitor (IFerrierd 
120011) . 



2.2. Mean-field model 

We follow the common paradigm where the magnetic field B 
and the velocity u are separated into mean part B and u plus 
fluctuations denoted by B' and u'. In this framework, the am- 
plification of the mean field is described via an averaged induc- 
tion equation that comprises an additional term £ subsuming 
the action of the turbulent flow, i.e., 



d,B = Vx[(u + tfQ*y)xB + £-77VxB] , 



(1) 



where rj is the microscopic magnetic diffusivity and the term 
including qQ. describes the induction due to the background 
shear. 

From the direct simulations described above, we compute 
the turbulent mean electromotive force £ = u' xB' by applying 
horizontal averages. From the obtained vertical profiles £(z, f) 
and u z (z, t) we successfully reconstruct the mean magnetic field 
B(z, f) via Equation (Q3 demonstrating the applicability of the 
chosen mean-field approach. 

As a simple closure for the mean induction equation, that 
via £ still depends on the small-scale flow u' and magnetic field 
B', we assume a standard parameterisation of the form 
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Fig. 1. Dynamo coefficients for model H4. Starting at t — 
1.3 Gyr the TF-method is applied to four consecutive time in- 
tervals of 25 Myr each. Quantities indicated by the ordinate la- 
bels are plotted in dark (q-rr, . . . ) or light (ar^, . . . ) colours, 
respectively. Shaded areas indicate lcr-fluctuations. 



While the diagonal elements of a produce the dynamo pro- 
cess itself, its antisymmetric off-diagonal part is responsible 
for the vertical pumping y z - \{a^R - a^). Similarly, the 
diagonal elements of f\ are interpreted as turbulent resistivity 
?7t = \(T\rr + fj<Mi)> while its antisymmetric off- diagonal com - 
ponents can lead to fixj-type dynamo effects (Ra dlerl[l969l) . 
determined by S z = \(fj R<p - fj^ R ). 



£, 



aajBj 



fiijEj k ,d k Bi 



i,jz{R,4>},k = z, 



(2) 



with tenso rial parameters a , and fig. W e apply the test-field 



3. Results 

We here compare six models (see Table [TJ with varying SN- 
rate cr and rotation frequency Q, keeping the shear parameter 
method of Schrinner et al ] (120051 12007) to invert this tensor q and the ratio cr 1 : cr 11 fixed. From the initial configuration 

turbulence builds up smoothly within about 50 Myr and the 
disk reaches a quasi-stationary state. The vertical stratification 
is now determined by the additional kinetic pressure from the 
SNe. 



equation for given test fields and their associated EMFs. In 



particular, we use the fields proposed by Brandenburg (2005) 



for the shearing bo x case (see also iBrandenburg et al] 
IGressel et al.ll2008l) . 
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Fig. 2. Profiles of the regular radial and azimuthal field for model H4 at various times. Results of the simulation (grey lines) are 
compared with the fields computed from &{z, t) via eq ([TJi (black lines). At t =s 0.85 Gyr, a field reversal with pronounced dipolar 
symmetry occurs. 
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Table 1. Overview of conducted models. The letters 'Q', 'H', 
and 'F' indicate quarter, half, and full SN-rate while numbers 
give the rotation rate. All models include shear with q = — 1. 



3.1. Dynamo parameters 

The most distinctive finding persistent in all our runs is that 
despite the peaked distribution of the SNe the velocity disper- 
sion increases with galactic height for the inner 2 kpc of our 
box. This is also reflected in the lower two panels of Fig. [TJ 
where we plot the a and f\ parameters obtained via the test- 
field method. In excellent a greement with the predictio n of 
second-order theory (SOCA, iRiidiger & Hollerbach"l2 004). the 
turbulent diffusivity closely follows the square of the turbu- 
lent velocity profile. From this comparison, we can infer a cor- 
relation time t of the turbulence of » 3 Myr. We measure rj t 
to be around 2 kpc km s~', with only a weak dependence on 
the applied supernova rate. Improving over previous studies, 
we determine the off-diagonals of f\ and find a positive value 
of 6 7 & 0.5 kpc km s" 1 . According to the dispersion relation 
in [Brandenburg (2005), a positive value for 6 Z will, in princi- 
ple, allow for a contribution of this term to the dynamo effect. 
Turbulent pumping is found to be directed inward and has an 
amplitude of |y, | « 5 km s" 1 . Most strikingly, in our simulations 
the inward pumping is balanced by a wind of approximately the 
same magnitude. Consequently, the effective vertical transport 
is drastically reduced, relieving the dynamo process of a major 
burden. 

3.2. Field structure and amplification 

In Fig. |2] we plot vertical profiles of the exponentially growing 
mean radial and azimuthal magnetic field for model H4. We 
find the field to be largely confined to the inner disk with a ra- 
dial pitch angle of -10°. Although this is somewhat smaller 
than the observed values of up to -30°, this value is con- 
siderably higher than would be expected from classical aQ- 
dynamos. While the predominant symmetry with respect to the 
midplane is found to be quadrupolar, this mode is interrupted 
by field reversals showing dipolar symmetry. This distinct be- 
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Fig. 3. Evolution of the turbulent and regular magnetic field for 
model H4. For (Br) and (B$) we show the results from the 
direct simulation (grey lines) together with the reconstruction 
from £(z, t) (black lines). 



haviour was successfully reproduced in a ID toy model, where 
the frequency of the periodic field reversal depended critically 
on the interplay of the diamagnetic pumping and the mean ver- 
tical velocity profile. 

To study the temporal evolution of the arising fields we in- 
troduce vertically integrated rms-values {Br) and ( B^ ). The 
exponential growth of the field is illustrated in Fig. [3] The e- 
folding time is on the order of 250 Myr and varies with the 
reversals. 

In our simulations, the turbulent component dominates over 
the regular by a factor of 2 - 3. In particular, we find values 
(B) : (B') of 0.52 (+0.02) for model Q4, 0.40 (+0.03) for 
model H4, and 0.31 (±0.01) for model F4. This trend is con- 
sistent with observations of st rong regula r fields in the inter- 
arm regions of spiral galaxies (Beck 2007). Furthermore, from 
IR-based star formation rates. lChyzyl (12007 ) observes a correla- 
tion log (_B reg : B tur ) oc -0.32 (+0.01) log SFR. From our values 
cited above we find a somewhat steeper slope of -0.38 (+0.01) 
with respect to the SN-rate. 

From the left panel of Fig.Hl where we plot the time evo- 
lution of the regular and fluctuating components of the vari- 
ous models, we see that even the absolute value of the mean 
field (B) increases with decreasing SN-rate (models F4, H4, 
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Fig. 4. Comparison of the time evolution of the regular and 
fluctuating magnetic field strength over kinetic energy for the 
various setups (cf. Tab. [TJ. Model F4* is identical to model 
F4 with Coriolis forces disabled and demonstrates that no field 
amplification is obtained in the case of Cartesian shear. 



and Q4). This is consistent with the trend in the turbulent dif- 
fusivity rj t , which has values of 2.3, 1.7, and 1.4kpckms _1 , 
respectively. 

For the range of parameters studied, we do not observe 
a significant dependence of the growth rate on the super- 
nova frequency <x but only on the rotation rate Q. For the 
models F1-F8, we find e-folding times of ~ 500, 140, 102, 
and 54Myr for the amplification of (£?')• With exception of 
mod el Fl, wh i ch di rectly corresponds to the parameters used 
in iKorpi et all (Il999h . we find exponentially growing regular 
fields (B) with e-folding times of 147, 102, and 52Myr for 
model F2, F4, and F8, respectively. For model Fl the regular 
field decays at * 500 Myr. The listed values have been obtained 
for a time frame of about 100 Myr after the turbulence reaches 
a steady state. Due to the different time base these values are 
not directly comparable to the long-term growth rate of model 
H4, where field reversals become important. The initial ampli- 
fication timescale for the models H4 and Q4 is 90 Myr. 

4. Discussion & Conclusions 

We have performed direct simulations of a supernova-driven 
galactic dynamo and find the rotation frequency to be the crit- 
ical parameter allowing the dynamo to operate. For our setup, 
it turns out that Q. > 25 km s" 1 kpc" 1 is necessa ry, which coin- 
cides with the prediction of lSchultz et alJ (fl994). Nevertheless, 
this value may still be dependent on the assumed gas density 
and the gravitational potential. 

While the e-folding time of the amplification mechanism 
scales with the rotation period, within the range of parame- 
ters studied, it only shows a minor dependence on the super- 
nova rate. With r e - 250 Myr, the overall growth time for our 
model H4 is compa rable to the value o btained for the cosmic 
ray driven dynamo (Han asz et al . 2006) and about four times 
larger than the expected growth time for the magneto-rotational 
instability. 

For similar models without shear dGressel et al.ll2008h . we 
found no amplification of the mean field. Further simulations 
will have to show whether the differential rotation is indeed 
essential for the SN-dynamo to work, or if it was only that 
the critical values for dynamo action were not achieved in the 
rigid rotating case given the assumed density profile and grav- 
itational potential. An estimation based on dynamo numbers 
leads to slightly subcritical values for mean-field dynamo ac- 
tion in the case of rigid rotation. 

Although the effects of numerical resolution have to be in- 
vestigated by further simulation runs, we are confident that the 



main features are indeed robust with respect to the numerical 
modelling. The good agreement of the direct simulations with 
the ID toy model, as well as the absence of field amplification 
in the case without Coriolis forces, provide further evidence 
in support of the classical picture of cyclonic turbulence. That 
is, the Coriolis force plays the central role in giving the ex- 
panding remnants a characteristic handedness. The mean ver- 
tical velocity, together with the turbulent pumping term, add a 
new timescale to the system, which is probably responsible for 
the fast growth. This was already observed in me an-field mod- 
els w ith a prescribed wind in the mean velocity dElstner et al] 
119951) . 

The relative strength of the regular field compared to the 
turbulent field is determined by the SN rate. The resulting 
slope of -0.38 (±0.01) is similar to the observed one derived 
from comparisons of polarised over total int ensity of radio syn- 
chrotron emission with SFRs dChvzy| |2007). The strong radio- 
FIR correlation also supports a strong dependence of the total 
magnetic field strength with the SFR. Therefore, the saturation 
level of the dynamo with respect to the SN-rate should be com- 
pared in the future. 

Since we do not observe evidence for catastrophic quench- 
ing, our results suggest that the combined action of the dia- 
magnetic pumping and the wind might cause an outward he- 
licity flux while retaining the mean field. This poses a natural 
solution to the quenching catastrophe and will be the subject of 
further investigations at higher magnetic Reynolds numbers. 
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